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CHAPTER I 


INTRODUCTION 

This work represents an attempt to characterize the essential dynamical 
characteristics of a simple single-spool turbojet engine through simulation • 
of low order system models on an Electronic Associates TR-.48 analog computer. 
The objective is to gain insight into the most important dynamical constraints 
of such an engine, which can be applied to control studies of more advanced 
engines. The approach here is to begin with an accurate model which is 
studied in [1] and reduce system complexity through various linearizations 
and approximations. References [2] and [3] have been used as guides. Some 
excellent ideas in [4'] have also been used, although this reference became 
available only after the present simulation was essentially in its final 
form, and thus certain suggested changes must be postponed for future, 
studies . 

This work consists of a derivation of a seventh order simplified 
simulation model, a derivation of an even simpler third order model, and 
simulation results from each. 

The control problem studied is one of getting from "Windmill” (zero 
fuel flow equilibrium) to "Design" (a high thrust equilibrium) while taking 
into account surge margin and turbine inlet temperature constraints. 

Several control schemes were investigated. 

As a matter of terminology, we wish to point out that in this work 
we. consider dynamical system models of the form 

X =» f (x,u) (1) 


1 



2 


with an associated design equilibrium point. 



such that 

f(xg, Ug) = 0 

We use normalized variables 


“i ■ 

1 


X. •= x./x„ 

1 1 E. 

X 


so that the normalized equivalent to (1) is 


( 2 ) 

(3) 

(4) 


X = F (x,G) (5) 

with an equilibrium, corresponding to (2) , occuring where components of 
Xg and Ug are all unity. It is also convenient to consider linearizations of 


the normalized system (5) through 

A /V 

X = 6x + Xg 

and the standard approximation 

A A 

6x = A6x + B6u 


■ where 


A A 


u = 6u + u. 


( 6 ) 


(7) 


A = ^ 

3x 

evaluated at (x_, u ). 

b b 


B = 3F 


3u 


( 8 ) 



CHAPTER II 


DERIVATION OF SIMULATION MODELS 


2.1 Introduction 

In this section we review the essential details in approximating the 
accurate simulation model of [1] by a seventh order dynamical system. The 
principal constraints lie in the limited nonlinear equipment available on 
the TR-48 analog computer. The main simplifications involved the extensive 
use of linear approximations, single stage compressor dynamics, a linear 
compressor map, assumptions of a choked nozzle condition, and certain ■ 
empirical relations based on design point equilibrium data available in [ 1] . 
Additional simplification was achieved by limiting the model to a condition 
of 20,000 ft. altitude at Mach number ,8. We shall term this model used in 
[1] the "Drone" system. 

2.2 The Drone System 

A complete description of the Drone system is presented in [ 1] , with 

certain details further elucidated in [2] and [4]. A rough simulation 

diagram is indicated in Figure 1. Throughout the development, P stands for 

2 

pressure in newtons/meter , T stands for temperature in degrees Kelvin, o) 
represents mass flow rate in kilograms /second, and p is density in kilograms/ 

3 

meter . N is rotor speed in revolutions /minute. Using a single stage rather 
than a four stage compressor model the main dynamical equations are: 


3 



4 


Compressor 


= 1 


c- 


Burner 


Turbine 


Rotor 


dt 

dP^ 

£, 

V 

''1 

= RY 

dt 


du)2 

- 

dt 



1 


= ^ c; 

dt 



= RY 

dt 


d(i), 

4 

= \ ( 

dt 


ddif" 

= i( 

dt 

T 

dPr 

= RT^ 


5 

dt 

''5 

^ = 



2 “ 3 ^ 


(T2<D2 - ^2' “3> 


4^3 ^ f • ~ 5^ 4 ^^4“B “'f 


dt \ TT 


fp 

NI ^ pB (iL^T^ 


( 1 ) 


( 2 ) 


(3) 


(4) 


(5) 


( 6 ) 


(7) 


( 8 ) 


Sc^“2'^2 


W 3 T 3 )' (9) 
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Auxiliary Relations 


P2- = RPJ2’ 


^^4 ^ ^ “4 ^^3 


(O3 = 0)^ 


Non-Linear Functions 
Turbine 

“3 ° ■ 

h ‘ * 2 <'' 4 ’'' 4 ’' 5 ’’'> • 

Nozzle 

“8 ^ 


( 10 ) 


( 11 ) 
, ( 12 ) 


(13) 


Compressor ' 

T 3 - f 4 (T 2 ',N,P^) (14) 

P3 = f^(P2',T2*, 0)3,102) ( 15 ) 

2.3 Assumptions 

We now list the assumptions which were made in order to reduce the required 
nonlinear function generating equipment. Where appropriate, we point out the 
corresponding simplifications suggested in [4]. As mentioned earlier, many of 
the suggested changes were not implemented simply because they were not avail- 
able by the time the simulations were run. One general comment is also in 
order with regard to linearizations. Rather than take partial derivatives of 
theoretical relations, the usual approach here was to take the best linear fits 
to the numerical data available for the five design points studied in [ 1 ]. 
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(a) = 280. 6°K. This eliminates .the need for- state equation 
(1) involving p^' and. is consistent with [4]. 

(b) This .ignores combustion lags and dead time. This assumption 
is also made in item of reference [4] . 

(c) For d)^ = f^(N,T^,P^ P^) related to turbine characteristics we have 

• -6 

used simply oj^ = 12.424 x 10 P^ based on numerical data fits. This is in' 
contrast to item of [4] which suggests a simplified turbine map to obtain 



(d) For d)Q = f„(Tc>F(-) have used co„ = 0k„A„P,./ i^tT = 22.835 x.lO ^ P^0 

ojoo ooo55 5 

which amounts to an assumption of choked noz 2 le condition. This is exactly 
the same as (15^ of [4] except that rather than linearize 1/*'*^ we finally just 
use a constant approximation. 

(e) For T^ = f^CT^jP^jP^jN) we simply use the empirical relation 

T^ = .87694T^. This corresponds to item ( jL?) of [4] which suggests linearizing 
the theoretical relation to obtain an approximation ' of the form T^ = A + BT^ + CP^. 

(f) For T 3 = f^(T 2 ',N,P^) we use T 3 = T^' + 12.914 x 10“V which is a 

form suggested by Zucrow [3], but differs with item of [4] which suggest a 

linearization of the theoretical relation to express = A + B P^. 

(g) For ?3 = f^CP^' ,W 3 ,d) 2 ) which is the compressor map relation we have 
used essentially the method suggested in item ® of [4] with the constant speed 
line function simply a straight line. This results in P^ = 5*10 ^P^ N - 2.94104 
X lO^u)^. 

The above assumptions yield a seventh order nonlinear system which 
will be shown later in this chapter. Additional assumptions made to obtain a 
third order simplified model are: 

(h) Replace equations (2) and (3) by their corresponding equilibrium 
relations. Thus, we neglect compressor dynamics and bleed losses, corresponding 


reproducibility of the 
ORIGINAL PAGE IS POOR 
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to items (I) and @ of . 

(i) Replace equation (6) by the corresponding equilibrium condition 

relates, to item, (s) of [4] which lets = 1.05 P^. 

« • 

(j) Replace equation (8) by the equilibrium condition a:^ = which 
results, with our above assumptions, in the condition P^ = (22. 835/12. 424)P^0. 
This assumption is not made in [4] , accounting for the fourth order system 
there and the third order system in this work (We use 0 < 0 < 1 as a 
parameter for nozzle area) . 


2.4 Thrust Calculation 

Here we follow Zucrow [3^ to compute net thrust due to internal flow. 
We have 


F . (Vj - V) + (Pg - P^) Ag 


(16) 


or 


F = F. (S- V 
3 


(17) 


where 


■ ^3 - V!i/3 ^8-^ (P - P ) 

g J 0)8 8 o 


(18) 


and for the choked nozzle condition 

Y 


Y 




(19) 


Which becomes for Ag = .01820 and y ~ 1.4 


F. = 0 P^ - .0182 p„0 


10 ' 


10 


( 20 ) 


With the Mach .8 at 20,000 ft. condition we get g = 1, V = 252.89 m/sec, 
pQ ” ^‘86 X 10^. Also we have by assumption (d) , tlig = 22.835 x 1O~^P^0. 
Using an efficiency onY) = .9464 to fit the design data of [1], and equation 
(17) we obtain 
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F = 0 {.01658 Pg - 848.12} 


( 21 ) 


as the approximate expression for- thrust.. 


2.5 Seventh Order Simulation Model 


State Equations; 


" 11-2745(10^) - i^) 

dt“ \ 


( 22 ) 


dW 2 (7.09(10^) - ^ 2 ’^ 


dt 

’1 


da .3 


P 3 - ^ - 

dt 

^4 



= 401. 

8 [’^4“f + 

4 

dt 


287p^ 


= 1 

[ OJ -f-w- - 

dt 


J z 


4.108(10^) m 

"3 


“f ^ ‘“a "3 


12.424 P, 


(23) 


(24) 


10 P 


B 


(26) 


10 ' 


dP 

dt 


5 .87694 P, 


,„ 6 „ (12.424P. - 22.835P.0) 

10 V 5 Pg 4 5 


(27) 


dN = 1.0616(10^) r ^4 .12.424 20.025 , 

dt IN ^ 287p^ ^ ...6 ^4 " ,.6 V-* 


B 10 


10 ' 


+ .8633 [280. 60)2 - T^Wg]} 


(28) 


ItEPRODUCIBILri’.*' . 
OBIGMAL PAaS 1 
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Other relations; 


='1164.15' = 1005.02 'yi4lig/Cpg =-36278' 

= 280.6 J = 1 Kg = 4.108 x 10® = 7.09 x lO'^ 

0 ). = i, P, = Rp T, A„ = .0182 R = 287 y = 1.4 

O 4 4 15 4 o 

“5 = 12.424 X 10“® P, 

4 

o)o = 22.835 X lO"S-0 

Tc = .87694T, 

5 4 

Tg = 280.6 + 12.914 X 10“®N^ 

Pg = 5P^N X 10“'^ - 2.94104 X 10^ 

F = 0 '{.Ol&SSP^ - 848.12} 


( 29 ) 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


2.6 Gains 

It was impossible to obtain all the gains determined by the factors: 

2 2 

1/yf^, /V^, 1/V^, I/V 3 , and l/l in the equations above; although some 

information given in [4] was useful. To get a range of dynamical characteristics, 
four gain systems were considered: 


System A: 

Each V = 

.1, each A = 

.1, I = 

.00305 

System B: 

Each V = 

.1, each A = 

. 1 , I = 

.0305 

System C: 

Each V = 

.01, each A = 

II 

H 

.0305 

System D: 

Each V = 

.01, each A = 

H 

11 

.305 


Note, in referring back to the model of the previous section, that System A 
■is obtained from System B by multiplying equation (28) by a factor of 10. 
System C is derived from B by multiplying equations (22)’ - (27) by 10 and (28) 


Tii' 
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by 1. ' System D Is obtained from System B by multiplying (22) - (27) by 10 
and (28) by 1/10. 

2.7 Normalized Seventh Order Simulation Model 

Here we use design point (2) of [1] for the equilibrium values. The 
associated normailzed equations then follow as in the introduction. Gain 
System B is used here. 

State Equations: 


dPl . . ■ 

IT = ^“2 


do), 

= 2188.27 (1 - Pp 


^^^3 = 8191.36 [1.05787P - P, - .05787 
dt P3 


.93586 ^4"f , 

i- 31.486(1)- + 21.435u)oT_ 

dt Pg f 3 3 


53.86 ^4 


dp „ r. 

= 37.78(U3 - 38.448 + .66849m^ 


dPp ■' 

dT = jr t^4 - V' 

D 


" P ^ 

4^ = — ^ {-A [3.12P, - 2.7361P.0] 
.305N ^ 


+ [.688013(L2 - 1.0715 T^W3l> 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 


Other Relations: 


A, A* A A A 



( 43 ) 
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y\ 

5 4 


A 



T 3 = .64212 + .35788N^ 

= 4.394N - 3.394(02 

F = 0 {I. 5486 P 2 - .5486} 

V = V, = = A = A, = .1 I = .0305 

14 5 14 

(this corresponds to Gain System B in 2.6) 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


2 . 8 Normalized Third Order Simulation Model 

Here we use assumption (h) , (i) , and (j) of Section 2.3 to eliminate the 
high frequency terms of the seventh order model. Note that the equilibrium 
solutions of both systems are identical. Here again we use Gain System B. 
State equations: 


dP 


dt 


^ (.93586P^/Pg + 31.486) + 21.435to,T, - 53.86 P^/p, 


3^3 


■4 


dp ^ y. ^ 

3 — 2 - 37.78(0,, - 38.448P, + .66849co^ 
dt 3 4 f 


dN 1.258 2,.^ : „2, 

d£ ‘ ~ir - “ 3 “ > 


(51) 

(52) 

(53) 


1. 
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Other Relations: 

All relations equations (43) to (50) in the previous seventh order 


model are valid, plus 

P2 = 1 “2 “ ‘"S 


^4 = ^ 5 ® 


(54) 


0)^ - 1.3009N - . 13982 (P^ v'^^2 ^43^533^2 _ ,0899P^N) 


(55) 


2.9 Equilibrium Conditions 

Equilibrium conditions for the normalized seventh order and third order 
models are the same, and independent of whether Gain Systems A, B, C, or D 
are used. The "Design" equilibrium occurs when all normalized state variables 
are unity, and this corresponds to the design point (2) of Mach .8 and 20000 
ft. in [1]. Other equilibrium points must be calculated using a successive 
approximation procedure to solve the nonlinear equations which result by 
setting all derivatives zero. In particular, we define "Windmill” as the 

? A 

equilibrium which occurs when fuel flow = 0 and nozzle area 0 =1. 

The algorithum used to calculate equilibrium conditions is as follows. 

A 

1. Set values of and 0 (the controls.) 

2. As initial estimates, set 

P 2 = 1, N == P^ = .5(0)^ + 1) 

3. Set in order 

= 1.3009N - .13982(P^ - 2 _ .0899P.N) 

4 H 

T^ = .64212 + .35788N^ 

A A A 

P^ = 4.394N - 3.394u2 

P^ = (i^ + 56. 515^3)757. 515 
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= (53.86P^^ - .93586P^a)p/(21.435u^T2 + 31.486o)^) 


A A 


T = P /6 
4 4''^B 


Pj = p,/0 


4. Set 


N = 2,^ : 

P 4 /Pb1*3 

and return to 3. until convergence is achieved. 

This algorithm gave five place convergence in twenty iterations and showed 
a very nearly linear operating line relation of approximately - 

A A 

P^ = 1.0263w^ - .0263 (56) 

or equivalently (using the linear compressor relation.?, = 4.394N - 3 . 394 a) 2 ) 


0 ). = . 99405N - . 00595 


(57) 


Some values of interest are tabulated in Table I for 0=1. 

2.10 Compressor Map and Surge Lines 

The equation used to approximate the compressor map of Figure 3. in [1] 
is 


where 


Pj/F^ = 5N/10^ - 2.94215 o)2v^2^^2 


02= T 2 / 288.3 = 280.6/288.3 = .97329 
62 = P2/10-1325 X 10^ 


Using these relations and the equilibrium values 

P_ = 28.076 X 10^ 

3E 


P2g = 7.09 X 10 


“SE = 3*24 

vr ^ ^ 


(58) 


(59) 


( 60 ) 



TABLE I 


SAMPLE EQUILIBRIUM VALUES 


r<, 


“£ 


N 

•^B 

T 

4 



0.00 

.53831 

,54589 

1.77504 

.30326 

.54783 

.53931 

0.25 

.72430 

.73097 

1.34002 

■ .54051 

.73269 

.72513 

0.50 

.83809 

.84284 

1.17141 

.71545 

.84407 

.83867 

0.75 

.92614 

.92861 

1.07041 

.86522 

.92926 

.92640 

1.00 

.99998 

.99997 

1.00004 

.99994 

.99998 

.99994 


.74876 (Windmill) 

.83334 

.89635 

.95073 

.99998 (Design) 
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we obtain 

Pg = 4 . 394 P 2 N - 3. 394*3 (61) 

This relation- is graphed in Figure 2. for P^ = 1 along with the operating line 
of equilibrium points calculated in the previous section. It would be natural 
to map a linear approximation of the surge line of Figure 3 in [1] on to Figure 

I 

2, except' that our Windmill eqiulibrium, unfortunately, does not match the 
image of the windmill eqiulibrium in [1] . This is due to the many approximations 
made in section 2.3. 

We have thus chosen to define the surge line as a line running parallel 
to and above the operating line. This is rather arbitrary, but preserves 
the qualitative nature of the problem of getting from Windmill to Design 
without transgressing on the surge characteristic. Indicated on Figure 2. is 
a surge line specified by the relation 

P 3 = 1.0263*3 + .24105 (62) 


2.11 Linear Normalized Systems 

In order to estimate the dynamical modes of the. various models, and to 
provide models for linear control studies, linear representations of the form 


dSx 

dt 


= A6x + B6u 


(63) 


were determined, where A = ^ and B = 9f . Using the seventh order normalized 

9x . 9u ' 

equations of section 2,7, with Gain System B and linearizing about the design 
equilibrium 

Xg = (1, 1, 1, 1, 1, 1, 1) 


we have 
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A 

“2 


» 

A 

‘’b 

51.52 

-51.52 

0 

0 

0 

0 

0 

0 

0 

-31967 

-8191.4 

0 

0 

21.436 

-106.72 

52.893 

0 

37.78 

-38.448 

0 

0 

0 

61.97 

0 

2.2558 

-3.5131 

11.488 

-i:2587 
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The corresponding Gain System B third order normalized linearizations 


are 


A = 

A, 

N 


"4 

-112.27 

-48.12 

2.8377 


52.924 

0 

-1.258 


N 

42.26 

47.444 

-4.096 


( 66 ) 


B = 



N 


32.4219 

.66849 

0 


(67) 


With this A matrix we were able to determine the A matrix of the other 
Gain Systems rather easily. We were also able to calculate the eigenvalues 
of these matrices with the aid of a digital computer program. These results 
are presented below. 

Gain System A. (V = .1, A = .1, I = .00305) to obtain A and B matrices, 

multiply bottom rows of all matrices by 10. 

7th order eigenvalues: -31852, -141.43, -60.278 

(t ^ .0425) -44.907 ± j329.24 

-8.5923 ± j21.913 

3rd order eigenvalues: -89.288, -45.638, -18.304 

(t .0546) 


Gain System B. (V= .1, A= .1, 1= .030^ A and B matrices exactly as shown 
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in equations (64)- (67). 

7th order eigenvalues: -31892, -88.732, -61.218 

(t ^ .296) -25.204, -3.3818 

-33.738 ± J333.68 

3rd order eigenvalues: -81.219, -32.309, -2.8386 

(t -352) . .... 


Gain System C . (V = .01, A = .1, I = .0305) to obtain A and B matrices 

multiply all rows except the bottom one by 10. 

7th order eigenvalues: -318962, -814.82, -618.54 

(t ^ .339) -310.44, -2.9529 

-325. 38± J3339.7 

3rd order eigenvalues: -807.68, -316.2, -2.9166 

(t 2; *3^3) 


Gain System D. (V- .01, A= .1, 1= .305) to obtain A and B matrices multiply 

all upper rows by 10 and the bottom one by 1/10. 

7th order eigenvalues: -318966, -806.7, -619.58 

(t ^ 3.423) -316.34, -.29218 

-324.17 ± j3340 

3rd order eigenvalues: -807.25, -315.57, -.2924 

(t 2; 3.420) 


^iPRODUCEBILErY OF i 
■^iRIGiNAL PAGE IS P(.'.' 



CHAPTER III 


SIMULATION AND CONTROL STUDIES OF THE SEVENTH ORDER SYSTEM 

3.1 Introduction 

This chapter presents the analog simulation study of the seventh order 
Drone engine dynamics. of section 2.5. The purpose of the simulation was to 
develop a low order realistic representation of a turbojet engine which 
would coincide with the data of reference [1], With this model, we were able 
to gain further insight into the .equations which govern engine operations, 
as well as investigate various fuel controls. Gain Systems B and C were 
considered. (Recall the discussion of section 2.6.) Four different fuel 
controls were examined, and the results appear in the form of trajectories 
plotted on the outline of an approximate compressor map and time response 
curves. These figures are located in Appendix A. 

3.2 Simulation 

The actual simulation was conducted at the University of Notre Dame's 
Analog Computer Laboratory, An Electronic Associates TR-48, and a few small 
TR— 20 and TR— 10 analog computers were used. Patching diagrams were obtained 
with the aid of ANSIR 3, a digital computer program which generates all 
pertinent analog simulation information, given the differential equations. 
Details of the simulation are given in Appendix C. 

A preliminary check of the seventh order system was conducted to verify 
its consistency with theoretical equilibrium values. The ratio of Windmill 
condition to Design point values of the state variables was calculated. The 
same ratio was also computed from experimental data. Table II, illustrates 
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these findings. Discrepancies between the two sets of ratios arise from the 
inaccuracies in the equipment used, especially the nonlinear function 
generators. 

TABLE II 

SEVENTH ORDER SYSTEM CHECK 


VARIABLE 

CALCULATED RATIO 

EXPERIMENTAL RATIO 

^2 

1.000 

1.000 

“2 

.5478 

.5752 

“3 

.5478 

.6000 


1.774 

1.655 

^4 

.5383 

.5854 

^5 

.5388 

.5237 

N 

.5459 

.5831 


3.3 Control Problem 

As stated earlier, the control problem is to schedule fuel in order 
to accelerate the engine operating state from Windmill to Design equilibrium 
conditions. However, maximum acceleration potential is limited by compressor 
surge or stall. As fuel is increased to accelerate the engine, pressure ratio, 
P 3 /P 2 , Increases and airflow, decreases. When stall is reached the pressure 
will drop and will cause a similar decrease in compressor efficiency. It 
is desirable to design controls which will cause the compressor to approach 
the stall line however, since these controls tend to give a faster engine 
response. Trajectories were plotted on the outline of an approximate compressor 
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map to provide an illustrative look at surge margin constraint. Since this 
thesis is only an introductory study of the Drone engine, the surge line was 
considered arbitrary. Although, the underlying thought in this investigation 
was to regard the stall line as minimal. 

Several coordinate systems were tried to observe the effect of the control 
on surge margin. These curves were determined using Gain System B and open 
loop control.' The first natural choice was to use the pressure ratio versus 
airflow as shown in Figure 3(a). Figure 3(b) follows from the assumption 
that compressor inlet airflow equals outlet airflow (oj^ = m^) . - Time response 
curves of indicate that it remains constant except for slight oscillations 
,at the outset of the transient. By assuming constant, we obtain the plots 
in Figures 3(c) and 3(d). We also chose to make a graph with and N as 
the coordinate axis in Figure 4. This was done so that surge margin studies 
of the third order model could be carried out using state variables. We chose 
versus and P^ versus N as our representative coordinate systems. 

3.4 Fuel Controls 

Four different fuel control concepts were investigated on the simulation 
with Gain System B and C. They consist of an open loop control, linear 
feedback design, and two nonlinear controls developed by empirical methods. 
There effects on engine performance and surge constraint are indicated by 
time response curves „and traj.ectory plots in Appendix A. Two equations 
for each control are given. The first equation being the control used for the 
d 3 niamic equations of 2.5 and the second pertaining to normalized equations 
of 2.7. The value of at the design point is one. 

Control 1 - Here we use the term open loop to refer to this control, 
it simply involves letting 

0)^ = .057353 kg/ sec (1) 

EEPRODUCIBILrUY OF 1- 

_ ^ Jt TT1 tC* t 


or 
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i^= 1 ( 2 ) 

This corresponds to the value of at design equilibrium. State variable time 
response graphs for System B are in Figures 5 and 6 and System C graphs 
in Figures 7 and 8, 

Control 2 - Control 2 fuel scheduling is analogous to Case B control 
of [1], It is a linear combination of compressor and turbine pressures given 
by 

= 2.02954 X 10~^ (P^ + .032047Pp kg/sec (3) 

= .99392 + .00608 P^ (4) 

This control provided another means of comparing our simulation results to 
the data of reference [1] . System B acc-eleration transients for this control 
are in Figures 9" and 10 while System C curves are in Figures 11 and 12. 

Control 3 - This control is the product of compressor airflow and 
rotor speed. It is defined as 

= 5.0867 X 10“^ (m^N) kg/sec , (5) 

A A A 

N (6) 


Figures 13 and 14 give time plots of Gain System B and Figures 15 and 16 are 
of System C. 

Control 4 - Figures 17, 18, 19, and 20 show engine dynamic time response 
for this nonlinear control given by 


d). = 6.982 X 10 ^^(p//p^) kg/sec 

X H ^ 


“C = P/. /P, 


“f ‘4 '^B 


(7) 

( 8 ) 


3.5 Results 

Relative time responses for each system and control are presented in 
Table III. The time constants were determined from the response of the rotor 


speed. 
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TABLE III 

TIME CONSTANTS (sec) 


FUEL CONTROL 


SYSTEM B 


SYSTEM C 


1 

2 

3 

4 


,196 

,226 

.500 

,426 


.266 

.432 

.633 

.533 


These agree fairly well with the linearized results of section 2.11. A 
comparison to the Case B acceleration transients of [1] to our corresponding 
control 2 and System C. System C was chosen since it most closely represents 
the engine described in [1] and [4]. Our time constants are approximately 
.432 sec for N and .456 sec for Whereas the results in [1] indicate 

response times of .5 sec for N and .565 sec for 

Analysis of these four controls on compressor surge indicate that 
Control 1 and 2 allow for approximately the same surge margin with 1 yielding 
a slightly higher compressor surge. Controls 3 and 4, however, would tolerate 
a stall line much less than the one required by 1 and 2 to insure efficient 
compressor operation. This is particulary evident in System B, Figure 21 where 
we were able to plot the trajectories of all four controls on the same graph 
without confusion. System C yields the same results as shown in Figures 22 
and 23. 

We note that although controls 3 and 4 tolerate lower stall lines, they also 
require more time to accelerate the engine from Windmill to Design. 


3.6 Comment 

Analog computer results seem to suggest that the seventh order model is 
a good representation of the Drone engine. The major difference between 
Gain Systems B and C is in the response time, with System B yielding a faster . 
acceleration time. However, it appears that System C most closely approximates 
the theoretical engine design in [1] necessary for further control studies. 



CHAPTER IV 


SIMULATION, AND- CONTROL STUDIES OF THE THIRD ORDER SYSTEM 

4.1 Introduction 

The third order Drone engine simulation results are contained in this 
chapter. The purpose in developing this model was to examine the possibility 
of achieving an even simpler representation of a turbojet engine and still 
retain the essential characteristics of engine dynamics. Gain Systems A, 

C, and D were used in this study. The four controls discussed in section 3.4 
•were also implemented on this simulation. Slight modifications were necessary 
in the control schemes to be consistent with the third order simplifications. 
Also the effect of a linear combination of Controls 1 and 4 was investigated 
on surge margin constraint. The normalized equations of 2.8 were simulated 
and details of the actual simulation are in Appendix D. 

4.2 Simulation Check 

A check of the third' order model where ratios of Windmill to Design 
point equilibrium values were calculated, reveal's that the third order system 
more closely approximates the theoretical values. The major discrepancy in 
0)^ is probably caused by- the dynamical equation characterizing its behavior. 
The airflow is given as an equality rather than a differential equation. 

Also, the numerous simplifications needed to reduce the seventh order model 
undoubtedly had a significant effect. The results of this check are shown in 
Table IV, 
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VARIABLE 


Ur 


N 


TABLE IV 

THIRD ORDER SYSTEM CHECK 


CALCULATED RATIO 

.5478 

1.774 

.5383 

.5459 


EXPERIMENTAL RATIO 

.4348 

1.852 

.5324 

.5472 


4. 3 Fuel Controls 

The four controls studied were basically the same as in 3.4 the only 
.change occuring in Control 2. With the simplifications necessary to derive 
the third order model. Control .2 reduces to 

% = ^4 (1) 


The third order normalized fuel controls each have equilibrium values 

of 1 and are given by 

Control 1 = 1 

r 

' " ( 2 ) 

Control 2 = F, 

f 4 


Control 3 

r 3 

? " 2 

Control 4 = P. /p 

r 4 i5 


(3) 

(4) 


4.4 Results 

Time response curves for Control 1 are shown in Figures 24 and 25; 

Control 2 in Figures 26 and 27; Control 3 in Figures 28 and 29; and Control 
4 in Figures 30 and 31. The first figure number for each control refers to 
System A response and the second to System C. Time Constants measured on these 
plots are given in Table V. 
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TABLE V 

TIME CONSTANTS (sec) 


EUEL, CONTROL 

SYSTEM A 

SYSTEM C 

1 

.049 

.299 

2 

.059 

.466 

3 

.119 

.799 

4 

.153 

.633 


Time constants determined from the linearized model of .0546 sec. for 
System A and .343 sec. for System C agree fairly well with the experimental 
time constants of ,049 sec and .299 sec. 

Surge margin studies of Systems C and D yield the same conclusions -as 
in the previous Chapter. Controls 1 and 2 would be preferable if the stall 
line was sufficiently distant from the operating line. As the surge line 
approaches the operating line Controls 3 and 4 would tend to give better 
engine efficiency. These plots are shown in Figures 32, 33, 34-, and 35. 

Gain System A, however, demonstrates the necessity for improvements to be 
made in the simulation model. While System A still perserves the same qualitive 
results, we see that Control 4 would allow for a slightly smaller stall line 
than Control 3 as shown in Figure 36. This is contrary to what we have observed 
from the other Gain Systems. The reason being the large gain on the rotor 

/Ni .N /% 

dynamics, equation causing = ui^N to respond faster, which in turn gives 
rise to a somewhat larger pressure surge at the outset of engine acceleration. 

4.5 Linear Combination Control 

From trajectory plots that were determined from the simulation, it was 
decided to test a linear combination of Controls 1 and 4. The Intuitive idea 
was to design a control which would be time optimal (Control 1) with a 
control which would tolerate the greatest surge margin. The form of the 


control is 
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= K^(l) + K2 (P4^/Pb) . (5) 

where and were varied between zero and one; and = 1- 

Representative trajectories shown in Figure 37 were obtained using Gain System 
A. This controller is pleasing in the sense that almost regardless- of the 
location of the surge line, a control of the form of equation (5) could be 
designed so that compressor stall would not occur. The lower limit on the 
surge line constraint and still have this control avoid engine stall is the 
case where = 0 and equation (5) reduces to equation (4) . We also might 
note that the time to accelerate from Windmill to Design increases as K 2 
increases. 

4.6 Comments 

It should be- pointed out that the equations which characterize rotor 
dynamics (42) for the seventh order model and (53) for the third order 
system) both involve a division by N. Due to the lack of available nonlinear 
equipment for the seventh order model, this was set equal to the constant 
value of N at Windmill condition. The third order model allowed for the 
division of to be carried out. The effect of the approximation made on the 
seventh order model is not really certain, although the results seem to 
indicate a valid assumption. 

The third order model results lead us to believe that our simulator 
is a good representation of the Drone engine. We were able to reduce the 
number of nonlinear elements such that the entire simulation could be conducted 
on the TR-48 and on TR-20 analog computer. This reduces the inaccuracies 
often times encountered in large simulations without significant effect on 
the results. This model will provide a quick, easy, and inexpensive method 
of developing and testing time optimal controls before they are attempted 
on a digital computer. 
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CHAPTER V 
CONCLUSION 

This thesis represents an introductory study into the Drone engine. 

As such, the results seem good. It is unfortunate that all suggestions of 
reference [4] could not be implemented in this study, since the simulation 
was well underway before this information became available. 

Suggested further research regarding the simulation is to pinpoint exact 
values of volumes, areas, and lengths which geometrically describe the engine 
under consideration. Also the rotor inertia should be determined since there 
was an inconsistency in [4]. Placement of surge line and constraints on 
turbine inlet temperatures should be defined in order to determine an efficient 
control. 

Control design studies to be conducted in the future involve developing 
a time optimal control. This design will center around Dynamic Programming 
solutions to the Hamilton Jacobi-equation and Fletcher-Reeves Conjugate Gradient 
method for solving the Hamilton-Jacobi Canonical equations. With these digital 
computer results, a control can be tested on the third order model and hope- 
fully provide a control for more advanced turbojet engines. 
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APPENDIX A 

SEVENTH ORDER SIMULATION ACCELERATION PLOTS 


The following pages illustrate the acceleration transients of the seventh 
order system as well as graphs relating to surge margin studies. 













FIGURE 6. Acceleration 
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FIGURE 12. Acceleration transients, System C, Control 2 
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FIGURE 13. Acceleration transients. System B, Control 3. 
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FIGURE 16. Acceleration transients. System C, Control 3 






FIGURE 18. Acceleration 
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FIGURE 20. Acceleration transients. System C, Control 4 
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FIGURE 22. vs N coordinate system’ plots of surge margin studies using System C, 
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APPENDIX B 

THIRD ORDER SIMULATION ACCELERATION PLOTS 

The third order simulation acceleration transients and surge margin 
studies are illustrated in the following pages. 







FIGURE 26. Acceleration transients. System A, Control 2. 
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FIGURE 32. Graphs of surge margin studies of all four controls 
using System A. 
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FIGURE 33. vs. N coordinate system plots of surge margin studies using System C. 
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DETAILS OF SEVENTH ORDER ANALOG SIMULATION MODEL, 
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Control 2 
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FIGURE 40. Analog diagrams of the four controls considered in this study 
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TABLE VI 
POT SETTINGS 


POT 

VALUE 

01 

.3758 

02 

.1879 

03 

.0405 

05 

.0810 

06 

.6667 

07 

.6212 

08 

.4288 

10 

.8847 

11- 

.8862 

13 

.9627 

15 

.8700 

16 

.2954 

17 

.3333 

21 

.4132 

22 

.1333 

23 

.2421 

25 

.1200 

27 

.5612 

28 

.2500 

29 

.1960 

30 

.1000 

31 

.2667 

32 

.1200 

35 

.3750 

36 

.5000 

40 

.2662 

42 

.2956 

43 

.1976 

46 

.9132 

47 

.6697 

48 

.1933 

50 

.4978 

51 

.1294 

52 

.6608 

53 

.1000 

55 

.2000 

57 

.1112 

58 

.1333 
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TABLE VII 
STATIC CHECK' 


AMPLIFIER 

01 

02 

03 

04 

05 

06 
07 
10 
11 
12 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 

- 24 
25 

27 

28 

29 

30 

31 

33 

34 

35 

36 

38 

39 

40 

41 

43 

44 

46 

47 
A23 
A9 
B3 
C5 
D5 


■VALUE 


.0000 

.8860 

.0405 

-.0810 

-.4288 

.0810 

-.4288 

.8847 

.9627 

-.7737 

.8700 

TIKE 

-.8700 
■ -.8847 
..8740 

-.9627 

.5546 

.7747 

.8937 

-.9439 

.9439 

.0707 

.4679 

-.1280 

.1280 

.5400 

.0270 

.7709 

-.0405 

.4754 

.8213 

.0000 

.0880 

-.0140 

-.8213 

-.0065 

-.7169 

.2000 

.7169 

' .0140 
-.7568 

-.6843 

-.5540 

-.1786 



TABLE VIII 


. P. 


B 




SCALE FACTORS FOR SEVENTH ORDER MODEL 


Variable 

p* 

r» 


Scale Factor 
80000 


u). 


80 




40 


Beta = 30 


300000 


150000 


N 


40000 


500 


^3S 


.015 


P^P 50 /Pg 


,05 


600000 




.1 




.08 


EEKaODXJCKILrrY OF • 
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APPENDIX D. 

DETAILS OF THIRD ORDER ANALOG SIMULATION MODEL. 
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TABLE IX 
POT SETTINGS 


POT 

VALUE 

01 

.3936 

02 

.6733 

03 

■•5000 

05 

.5000 

06 

.^806 

07 

•.8356 

08 

.1250 

10 

.31^6 

11 

TIME 

15 

.1295 

16 

.2946 

18 

.6667 

21 

.1500 

22 

• .1155 

30 

.3211 

31 

.7158 

32 

.6667 

35 

.4722 

36 

.1180 


.5359 
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TABLE X 
STATIC CHECK 


AMPLIFIER 


00 

01 

02 

03 

04 

05 

06 
07 
12 
13 
16 
18 

19 

20 
21 

24 

25 
28 
29 

33 

34 
37 
42 
B3 


■^EHiODUCffiiLrry' of t 
mamAL page is pon- 


VALUE 


-,5000 
-.5000 
.5000 
.5000 
-.5000 
■ .1480 
.5000 

TIME 

.6667 

-.6667 

.5000 

.5000 

-.5000 

-.0500 

".3332 

.2500 

-.2500 

.5000 

-.5000 

-.3332 

.2500 

-.6666 

-.2500 

.2500 



« ^ <Q. < S < »3 < H < «3 


TABLE XI 

SCALE FACTORS FOR THIRD ORDER MODEL 


Variable 


Scale Factor 

.5 


■ .5 


.5 


.5 




.5 

.5 


1.5 


Beta = 3 
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■ APPENDIX E. 

SUMMARY OF' REFERENCE [4] 

A private comniunication from B. Lehtinen and K. Seldner dated June,,. 1975 
SIMPLIFIED DRONE SIMULATION w/VARIABLE NOZZLE 

i SEA LEVEL STATIC CONDITIONS (recommended conditions) 

• ■ P^ = 10. 13’ N/cm^ 

= 288.3 degK 
,®2 = = 1 

ii PARAMETER VALUES 

n = 86% (constant) 
c 

Ti, = 100% (constant) 

4 

= 85% (constant) 

5 

V^ = .006354 (.2244 ft^) 

- ■ V^ = .001099 (.03880 ft^)‘ 

I = .0305 N-m-sec^ (.0225 Ib-ft-sec^) 

h = 4.256 X 10^ j/kg 
■ c ■ 

C = .25 (4187.57) = 1046.9 
? 

Y =1.4 

iii OPERATING POINT DATA TO USE FOR LINEARIZATIONS 

a) Use Data given in Table I, column 1 of TMX-2537 [1]. 

Notes: 

a) Use SLS condition and un-corrected variables. 



b) • la the following write-up, "Eq D-" refers to a TMX-2537 equation 
number and "Eq J-" refers to a TND-6610 [2] equation number. 



L(h/cAR. i^ifr c»/*ec.-,r5 





FIGURE 43. Compressor Simulation 


© All dynamics in compressor are neglected. 

© f is one "representative" speed line obtained from Fig. D-3, the 


overall compressor map. 




Sp£iK.D LjfJS 


FIGURE 44. Compressor Map 


Fit the map of Fig. D-3 by linearly 
shifting the one "representative" 
speed line with N. "Bias" is adjusted 
to null out shifts when operating at 
the "representative" speed. 

Generate a linear approximation to 

the stall line, as shown (P ). 

m,.- . j T -^STALL 

Thxs xs used only as 

a stall indicator - 


Assume constant efficiency (n = constant) . Obtain T as a function of 


P^' using eq J-7 & J-8. 


T - T 
^3 2 


(P3/P2)' 


Linearize this relationship as shown. 

Obtain P^ from P^ assuming 5% pressure 
dropfP^ = 1-05 P^). 



FIGURE 45. Relation 
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FIGURE 46. Combustor Simulation 



Assume W 2 = W 3 (no bleed losses). 

Assume 0)3 = 104 (neglect momentum effects) . 

Ignore combustion lags and dead times in combustor 

• • 

Eg J-25 then becomes . 

( 9 ) Assume constant efficiency = constant) . 

Use the energy eq (eq J-24) . 

(PT^) = ^ [<^3 T3 + (m^ - m5)T^ + 

Include continuity (eq J-26) and state equations. 


>A - If 



FIGURE 47 . Turbine Simulation 
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li) Continuity eq (eq D-B14) 


*^^5 R T_ >• • > 

— = - 5 (0.5 - 0.8) 


12) Energy eqs 


T - T = 

^4 5 C 


.h 
A a 


(Eq D-B13) 


P. 2/7 

Ah. = C T. [1 - (:^) 1 (Eq D-B15) 

1 p 4 


Assume constant turbine efficiency; nS = constant 


Ah « R_ Ah. 
a 5 1 


Combining , we obtain 


T, - T_ P. 2/7 

Af — ^ - n,[l - (/) 1 

^4 4 


Linearize this temperature-pressure relationship, as was done for 
the compressor (item @ ) . 

Turbine map (f^) 
use eq D-Bll 




/T 


4 


UJ.^v.4 1 


Let be a representative corrected speed line from Fig. D-21. 


Assume a linear shift in this 
speed line with N/V^^ . Use a 
L/i>iiT VjiLxJi limiter on the output of the 

multiplier as shown on the analog 
diagram. 





% 


FIGURE 48. Turbine Map 
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about operating point. 


FIGURE 49. Nozzle 



Assume choked flow in Nozzle (0 ^ /P^ ^ , 530) 







Eq (D-B21) then becomes 


J5/ 


M 30 

dt ^ IT 


Jc 

_E. 

I 


“5 



+ - I 5 ) 



Nonlinear Complement is 14 multipliers and 2 DEC'S. 

To reduce multiplier complement to 6, linearize the multiplier 
functions denoted by in the diagrams; 
for example, on page 3; 


pT, 


I 


( UT, - 


^■^4 > 

o 


^42 

0 


P T, 
0 4 


(T, - T 


'‘o) 


4 

0 


— T 
T 4 

o 



Note that throughout the simulation diagrams, no attention has been 
paid to scaling or sign conventions . 
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